Berry phase in atom optics 
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We consider the scattering of an atom by a sequence of two near-resonant standing light waves 
each formed by two running waves with slightly different wave vectors. Due to opposite detunings of 
the two standing waves and within the rotating wave approximation, the adiabatic approximation 
applied to the atomic center-of-mass motion and a smooth turn-on and -off of the interaction, the 
dynamical phase cancels out and the final state of the atom differs from the initial one only by 
the sum of the two Berry phases accumulated in the two interaction regions. This phase depends 
on the position of the atom in a way such that the wave packet emerging from the scattering 
region will focus, which constitutes a novel method to observe the Berry phase without resorting to 
interferometric methods. 



I. INTRODUCTION 

The geometric phase [ll-Q manifests itself in many dif- 
ferent phenomena of physics ranging from the polariza- 
tion change in the propagation of light in fibers (4, via 
the precession of a neutron in a magnetic field t o 
the quantum dynamics of dark states in an atom [lO!]. 
The geometric phase has also been used in topological 
quantum com put ing (Tl| as realized, for example, with 
trapped ions [12] ■ In the present paper we propose a 
scheme to observe the geometric phase in the context of 
atom optics. 



A. Brief review of geometric phases 

The concept of the geometric phase arises in the con- 
text of a Hamiltonian which depends on a parameter 
which is slowly varying in time. When this variation is 
cyclic, that is the Hamiltonian returns to its initial form, 
the instantaneous eigenstate will not necessarily regain 
its original value, but will pick up a phase. This phe- 
nomenon has been verified in experiments with polarized 
light, radio waves, molecules, and many other systems. 
The most prominent example is the Aharonov-Bohm ef- 
fect [l^, which was observed in 1959 and interpreted [ij 
in terms of the geometric phase. Moreover, many familiar 
problems, such as the Foucault pendulum, or the motion 
of a charged particle in a strong magnetic field, usually 
not associated with the Berry phase, may be explained 
elegantly [2, ^] in terms of it. 

Since the landmark paper [l| on the geometric phase, 
extensive research, both theoretical and experimental, 
has been pursued on quantum holonomy 0, 3, a diabatic 
[l,i,[ll-[iland non-adiabatic HSIHI, cyclic [22, p"] and 
non-cyclic]^, [13] , Abelian |25| and non- Abelian 26., 23} , 
as well as off-diagonal 0, geometric phases. More- 
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over, geometric phase effects in the coherent excitation 
of a two-level atom have been identified [29-31]. Since 
geometric phases are rather insensitive to a particular 
kind of noise |32j , the y a re useful in the construction 
of robust quantum gates [ll|,i53 36]. Although proposals 
have already been given for the observation of the geo- 
metric phase in atom interferometry [lol . Ist} , so far only 
the dependence on the internal atomic degrees of freedom 
was investigated. In the present paper we extend this ap- 
proach by taking into account external atomic degrees of 
freedom, that is the center-of-mass-motion of the atom. 



B. Our approach 

We consider the scattering of a two-level atom from 
a near-resonant standing light wave formed by two lin- 
ear polarized running waves with identical electric field 
amplitudes and frequencies. The propagation direction of 
the two waves is slightly different and the electromagnetic 
field is detuned with respect to the resonance frequency 
of the atom. Within the Raman-Nath approximation 
[ssl Isgj on the atomic center-of-mass motion, adiabatic 
turn-on and -off of the interaction and with the rotating 
wave approximation, we obtain a condition for the can- 
cellation of the dynamical phase and show that the scat- 
tering process is determined solely by the Berry phase 
depending on the internal and external atomic degrees of 
freedom. The key observation in establishing this con- 
dition is the fact that the dynamical phase is antisym- 
metric in the detuning, whereas the geometric phase is 
symmetric. As a result, a sequence of two such scatter- 
ing arrangements which differ in the sign of their detun- 
ings eliminates the dynamical phase and leads to the sum 
the two corresponding geometric phases. To analyze the 
geometric phase we use the approach [40| based on the 
adiabatic eigenstates, that is the dressed state picture. 

Since the geometric phase imprinted onto the internal 
state is position-dependent, we propose a scheme to ob- 
serve the geometric phase based on the narrowing of the 
atomic wave packet. This application of the Berry phase 



might be useful in the realm of atom lithography |41j . 



C. Relation to earlier work 

It is for three reasons that our approach is different 
from earlier work on the Berry phase arising in the inter- 
nal dynamics of two-level atoms driven by laser fields: (i) 
in our scheme we compensate the dynamical phase; (ii) 
the geometric phase acquired by the internal states is im- 
printed onto the center-of-mass motion of the atoms, and 
(iii) our setup does not require a traditional interference 
arrangement. 

In a landmark experiment the dynamical phase of a 
neutron precessing in a magnetic field has been compen- 
sated by an additional 7r-pulse. In our scheme this 
cancellation of the dynamical phase is achieved by chang- 
ing the sign of the detuning as the atom interacts with 
the first and then with the second standing light field. 
Moreover, in the experiment described in 10] the Berry 
phase is observed in the internal states only and is read 
out by interferometry of these states. 

We extend these ideas to atom optics where the center- 
of-mass motion is treated quantum mechanically. Here 
we take advantage of the entanglement between the 
atomic states and the center-of-mass motion, which al- 
lows us to read out the information about the geometric 
phase using the dynamics and the self-interference of the 
wave packet. 
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Figure 1: Scattering of the wave packet '5 = '5 (a;) of a two- 
level atom by a standing electromagnetic field formed by two 
propagating waves of wave vectors ki and k2. The latter 
assume an angle a and n ~ a with respect to the i-axis. The 
field envelope along the y-axis translates, according to the 
relation y = vt, into the time-dependent function / = f{t) as 
the atom propagates through the field with velocity v — vey. 
The frequency uj of the field is detuned from the frequency of 
the atomic transition between the ground and excited states 
\g) and |e) by A. 



D. Outline of the article 

Our article is organized as follows. In Sec. II we formu- 
late the problem addressed in the present paper, define 
our model and evaluate the Hamiltonian describing the 
interaction of a two-level atom with an appropriately de- 
signed standing wave. Next we connect in Sec. Ill our 
model with the one of Ref. [Ij and recall the expressions 
for the dynamical and geometric phases. Here we em- 
phasize the connection between the dressed and atomic 
states. Since the geometric phase is determined by the 
path in parameter space, we construct in Sec. IV the 
circuit determined by the envelope of the electric field. 
Section V is devoted to the derivation of explicit expres- 
sions for the geometric and dynamical phases. In par- 
ticular, we consider the weak field limit where the Rabi 
frequency is much smaller than the detuning. As an ex- 
ample, we evaluate the geometric and dynamical phases 
for the Eckart envelope in Section VI. Sections VII and 
VIII are dedicated to the discussions of the cancellation 
of the dynamical phase and the read-out of the geomet- 
ric phase with the help of the center-of-mass motion. We 
summarize our main results in Section IX. 

In order to keep the paper self-contained we have in- 
cluded detailed calculations in several appendices. For 
example, in Appendix|X]we apply the WKB-method and 



perturbation theory to rederive the dynamical and ge- 
ometric phase for a two-level atom. Moreover, in Ap- 
pendix [B] we calculate the integrals determining the ge- 
ometric and dynamical phases for a special form of the 
field envelope, which smoothly switches on and off. In 
this case the path in parameter space circles many times 
around the origin. Finally, in Appendix [C] we evaluate 
the flux through these infinitely many windings. 



II. FORMULATION OF THE PROBLEM 

In the present section we formulate the problem of a 
two- level atom scattering off two running electromagnetic 
waves with almost opposite wave vectors. For this pur- 
pose we first establish the relevant Hamiltonian and then 
evaluate the matrix element of the interaction. We show 
that under appropriate conditions this quantity factorizes 
into a product of three terms, which allows us to imprint 
the center-of-mass motion onto a geometric phase. 



A. Set-up 

We consider the scattering of a two-level atom off a 
near-resonant standing light field created by two trav- 
eling waves of wave vectors ki = (fc cos a, k sin a) and 
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k2 = (— fccosQf, fcsina). Both fields form an angle a rel- 
ative to the X-axis of a Cartesian coordinate system as 
shown in Fig. [TJ The two running waves have identi- 
cal wave number, that is |ki| — |k2| = fc, but propagate 
against each other since (ki)^ — — {^2)^- As a result the 
electric field reads 

E(i, r) = Eo(r) [sin(k2r - cot) - sin(kir - cut)] , (1) 

where Eq = Eo(r) describes the position-dependent real- 
valued amplitude of the wave with frequency w. 

The time evolution of the state vector | ^E*) of the two- 
level atom interacting with the electromagnetic field E 
follows from the Schrodinger equation [H, 



where the Hamiltonian 



i/a-d-E I* 



2M 



(2) 



(3) 



is the kinetic energy operator of the center-of-mass mo- 
tion of the atom of mass M. Here, Ha denotes the 
Hamiltonian of the free two-level atom with energy eigen- 
states |e) and \g) and the corresponding energy eigenval- 
ues Ef, = hwe and Eg = hujg, shown in Fig. [T] The 
field frequency uj is assumed to be detuned from the fre- 
quency of the atomic transition between \g) and |e) by 
A = We — cjg — The interaction between the atom 
and the electromagnetic wave is described by the dipole 
moment operator d. Moreover, the carrot on E indicates 
the operator nature due to the quantum mechanical de- 
scription of the motion of the atom. 

The mean value of the velocity v of the atom in the 
direction of the y-axis is large and remains almost con- 
stant during the scattering process. For this reason we 
consider this motion classically, which allows us to set 
y = vt. 

In contrast, the motion along the x-axis is described 
quantum mechanicall y. M oreover, we make the Raman- 
Nath approximation [SSl . Is^ , that is we neglect the ki- 
netic energy operator Hem in the Schrodinger equation 
(12). Hence, the displacement of the atom along the x-axis 
caused by the interaction with the standing light wave is 
assumed to be small compared to the corresponding wave 
length. Since Hem is omitted in Eq. ([2|), the coordinate 
X is considered to be a parameter. Moreover, we neglect 
spontaneous emission due to the small interaction time 
T and the non- vanishing detuning A. 

As a result the Schrodinger equation ([2]) reduces to 



ih-\^)^[Ha-d--Ej\^). (4) 
In order to solve this equation we make the ansatz 

(5) 



where the amplitudes Ag and are functions of time t, 
but depend on the position vector r as a parameter. It 
is for this reason that we have dropped in Eq. ^ the 
carrot on E. 

Substituting the ansatz Eq. ([5]) into the approximate 
Schrodinger equation we arrive at 



dt \A 



with the Hamiltonian 



^ 1 / ftA V* 
2\ V ~hA 



(6) 



(7) 



containing the complex- valued coupling matrix element 

y(t;r) = -2e-''"V-E(i,r). (8) 

Here we have introduced the dipole matrix element p = 
(.g|d|e), which can be considered as real- valued in the 
case of the two-level atom. 

B. Matrix element of interaction 

The remaining task is to derive an explicit expression 
for the matrix element V defined by Eq. ([8]) in the pres- 
ence of the two running waves. For this purpose we repre- 
sent the sine functions in Eq. ([T]) as a sum of exponentials 

(9) 

and evaluate V neglecting terms oscillating with 2w 
which yields 



Vit;r) ~ -ipEo(r)(e 



— 'ik2r _ g — 2kir\ 



(10) 



Next we substitute the explicit form of the wave vectors 
ki and k2 into Eq. (ITU|) and find 



V{t;r)^2pEo{x,y)e 



— iky sin a 



sin(fcxcosQ!). (11) 



At this point we make use of the fact that the motion 
of the atom along the y-axis is treated classically and we 
can replace the y-coordinate hy y = vt. Moreover, for 
the sake of simplicity we assume that Eq is independent 
of X, resulting in Eo(x,2/) = Eo{vt) = £of{t), where 
/ — f{t) denotes the envelope function along the y-axis, 
as indicated in Fig. [T] The coupling matrix element V 
given by dTTI) takes then the form 



where 



V{t;r) = hVt{x)f{t)e- 



uja = kv sin a 



(12) 



(13) 



and 



il{x) = —p£oam{kxcoaa) = fio sin(fca;cosa) (14) 
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are the velocity-dependent Doppler and position- 
dependent Rabi frequencies, respectively. 

Hence, the coupling matrix element V consists of the 
product of three terms: (i) the position-dependent cou- 
pling energy hfl{x), (ii) the envelope function / — f{t), 
and (iii) the time-dependent phase factor exp{—iujat) due 
to the motion of the atom through the field. 



III. DYNAMICAL AND GEOMETRIC PHASES 

In the present section we connect the Hamiltonian Eq. 
([7]) together with the explicit expression for the matrix 
element V, Eq. (|12l) . to the Hamiltonian used in Ref. [l| 
to derive the geometric phase. This approach allows us to 
take advantage of the results obtained in Ref. pj] . More- 
over, we establish the connection between the dressed 
and the atomic states. 



A. Connection to Berry's approach 

For this purpose we make the identifications 

X = ReV, Y = lmV, Z = hA (15) 
and the Hamiltonian ([7]) takes the form 

^W-Kx + ,;y ""-T), (16) 

where the real- valued parameters X, Y and Z form the 
cartesian coordinates of the vector TZ. Here we have 
chosen the calligraphic letter rather than the normal one 
in order to bring out the fact that the vector IZ is not in 
three-dimensional position space but in parameter space. 
This Hamiltonian has the eigenvalues e*-^^ = ±e with 



^VX^ + Y^ + Z^ 



1 



and the corresponding eigenstates 



1 



f z±n 



(17) 



(18) 



where TZ = \1Z\. 

Now we return to the determination of the adiabaticity 
criterion. The condition for an adiabatic turn-on and -off 
of the interaction is that at any time t the rate of change 
of the light-field amplitude is much smaller than the 
spacing between the time-dependent quasi-energy levels. 
Hence, the parameters of the system should satisfy the 
following inequality 



V~dt 



(19) 



The expression Eq. (I12p for the coupling matrix element 
V allows us to cast Eq. (fT9|) into the form 



lc{f 
fdt 



« + p\^p{t) 



(20) 



For a function f{t) which decreases exponentially at large 
t, that is f{t — >■ ±00) oc exp(— jij/r), we arrive at the 
adiabaticity criterion 



|A|r » v/1 + (war)2 



(21) 



In the adiabatic limit, we remain in the adiabatic states 
1^''^='=^) which accumulate the phases in the course of time. 
When the time dependence of IZ is such that 7?.(i) re- 
turns at time t = T to its initial value 1Z{—T), the two 
instantaneous eigenstates |\E'*^*^) acquire [1] the dynami- 
cal and geometric phases 



and 



1 

i / em 



n-ds 



(22) 



(23) 



So 



where Sq denotes the surface of integration, i.e. the sur- 
face determined by a closed circuit forming during the 
one cycle of parameter change from t = —T to t = T. As 
a result the eigenstates after a cyclic change read 

|*(±)(r)) =exp(T^¥'^?)exp(T^</?s)|^'(±)(-T)). (24) 

The geometric phase is solely determined by the flux of 
the effective field 1^/271^ through the area enclosed by 
the parameter Ti. — TZ(t) during one period of the pa- 
rameter change, i.e. during the time interval 2T. 



B. Connection between dressed and atomic states 

The dynamical as well as the geometric phase are for- 
mulated in terms of the dressed states Eq. (|18l) which 
arise due to the atom-field interaction. However, in or- 
der to observe the geometric phase in an experiment, 
the atom should be prepared in a well-defined free-atom 
internal state. For this reason we need to connect the 
dressed states with the atomic ones, that is with 

\g) and |e). 

Before the atom enters the light field, that is at the 
time t ^ —T, the interaction V vanishes, leading to X ^ 
Y — and TZ = \Z\ ^ ^-l^l- Here we have assumed 
for the sake of simplicity that the envelope function is a 
mesa function with a sharp turn-on at — T and a sharp 
turn-off at T. Needless to say this assumption is not 
necessary and we will consider later a smooth envelope 
function. Due to appearance of the absolute value |A| 
of the detuning it is useful to consider the two cases of 
A < and A > 0. Indeed, for A < we find 



|vl/(+)(_r)) = and|vI/(-)(-T)) =-|e). 



(25) 



whereas for A > we arrive at 

|*(+)(-r)) = |e) and \^^-\-T)) = \g). (26) 
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Hence, in the case of A < the ground and excited states 
fohow 1^''^"'")) and l^"*^^-*), respectively, whereas, for A > 
the ground and excited states fohow and 

By using Eqs. ([5]), ([25]), and we obtain the 

total phase 



ipg = 13 signA + 7, 



(27) 



of the ground state acquired during the interaction time 
2T, where 



and 



P = cpD~\A\T 



7 = (p_B signA 



(28) 



(29) 



are the total dynamical and Berry phases, respectively. 

The same procedure results in the total phase of the 
excited state 



-^g-nei-A), 



(30) 



with G being the Heaviside function. 

Note, that we have neglected the phase contributions 
proportional to ujg and We arising from Eq. ([S]), since they 
are independent of the coordinate of the atom. Indeed, 
we shall show that the Berry phase as well as the dy- 
namical phase depend appropriately on the position and 
therefore can be detected by a narrowing of the atomic 
wave packet. 



IV. CIRCUIT IN PARAMETER SPACE 

In Sec. II we have derived the explicit expression Eq. 
(IT^ for the interaction matrix element V of the atom- 
field coupling. We are now in the position to discuss 
the path TZ. = 7?.(i) in parameter space traversed in the 
course of time. 

Since according to Eq. ([15]) the Z-component of TZ is 
given by the constant detuning A, the circuit lies parallel 
to the Xy-plane and is given by 



and 



X{t;x) = h\n{x)\f{t) cos{uJat) 



Y{t;x) = -h\n{x)\f{t)sm{ujct). 



(31) 



(32) 



Due to the position-dependence of ft the circuit depends 
on a: as a parameter. 

This curve is most conveniently described by the polar 
coordinates 



leading to 



p = h\n\f{t), (j} = io^t 



p^p(c^)^h\n\f{c^/u^). 



(33) 



(34) 




-r/2 



Til 





Y 































Figure 2: Translation of the electromagnetic field envelope 
/ = f{t) shown on the top into a circuit in parameter space, 
depicted at the bottom, illustrated here for the Eckart enve- 
lope Eq. (|35p . This function leads to an infinite amount of 
windings around the origin of parameter space. 



A model of relevant for an experiment relies on a 
smooth envelope, for instance, by the Eckart envelope 



fE{t) 



1 



cosh(t/r) 



(35) 



Although laser beams are usually modeled by a Gaus- 
sian profile, results obtained with the Eckart envelope 
are believed to be similar to those with a Gaussian 
one. Moreover, in the case of the Eckart envelope the 
Schrodinger equation ([6]) for the probability amplitudes 
to be in the ground and the excited state has an exact 
solution, which can be used to compare the results ob- 
tained within the Berry approach. 

We emphasize that in the case of the Eckert envelope 
the amplitude p given by Eq. (j33p vanishes for t = ±oo 
and consequently the circuit starts from and terminates 
at the origin. In the course of time the amplitude p first 
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increases and then decreases again. At the same time the 
angle increases monotonously. As a result the curve 
7?. = 7?.(t) defining the circuit circumvents the origin 
infinitely many times before it returns to it, as shown 
in Fig. [2l However, due to the increase and decrease 
of p not all windings will be visible. The number of the 
prominent loops is determined by the value of ujaT. 



V. EXPLICIT EXPRESSIONS FOR PHASES 

In the preceding section we have discussed the form of 
the circuit in parameter space dictated by the longitu- 
dinal mode function of the electromagnetic field. In the 
present section we calculate the resulting geometric and 
dynamical phases and analyze the weak-field limit. 



A. General case 

We start by obtaining explicit formulae for the phases 
due to an arbitrary but smooth envelope function. In Ap- 
pendixl^we rederive these expressions within the WKB- 
approach. 



Geometric phase 



According to Eq. (f23| we have to calculate the flux of 
the effective field 7?./ (272.^) through the area enclosed by 
the circuit in parameter space. Since the vector normal 
to this surface is in the opposite direction to the Z-axis 
of the Cartesian coordinate system we find that IZ- dS — 
-hAdS and Eq. ([23|) with the definition Eq. ([29]) reduces 
to 
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4r.|A| 



dX / dY 



1 



(X2 + Y^ + (ftA)2)3/2 



(36) 



In terms of the polar coordinates p and (p defined by Eq. 
([33]) this integral takes the form 



4>f p{4>) 



-h\A\ 



pdp 



(p2 + ft2^2)3/2 



(37) 



Here the upper limit p = p{4)) of the radial integral de- 
pends on the circuit in parameter space parameterized 
by the angle (p. The integration over runs between the 
initial 4>i = 4){—T) and final 0/ = 4>{T) angles of the 
curve. Their precise form is dictated by the shape of the 
circuit. 

The integration over p can be performed and with the 
help of Eq. and the new integration variable t = 

(b/uja we arrive at 



UJa_ 

2 



dt 



V|AP + |f7(x)P/2(t) 



(38) 



UJa_ 

2 



dt 



where 



(39) 



(40) 



is the dimensionless Rabi frequency. The position de- 
pendence of the geometric phase arises from the position 
dependence of a = a{x). 



2. Dynamical phase 

Next we turn to the dynamical phase. For this purpose 
we substitute the expression for V , Eq. (IT2t . into the one 
for the quasi-energy e, Eq. ()17|) . and find 



(41) 



Together with Eqs. and (gHl) for the total dynam- 
ical phase, we arrive at 



1 



(42) 



or 



/3=^y di(yTT^w)-i), (43) 

-T 

where we have recalled the definition Eq. (HU)) of the 
dimensionless Rabi frequency a. 

When we compare the expressions Eqs. (p8| and p2)) 
for the geometric and the dynamical phases we find the 
identity 



7 = 



dp 

d\K\ 



(44) 



which is reminiscent of the Kramers-Kronig relations. 
However, the WKB-analysis presented in Appendix A 
shows that Eq. is merely a consequence of a Taylor 
expansion in powers of Wa/IAj. It is interesting to note 
that in the formalism employed in the preceding sections 
this fact is hidden. 



B. Weak-field limit 

Finally we consider these phases in the weak-field limit, 
that is for a ^ 1, when Eq. (|43|) for the dynamical phase 
reduces to 



1 

/3 - i \A\a' J dtfit), 



(45) 
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while the geometric phase given by Eq. (p9)) reads 



B. Dynamical phase 



7 = -^t^afl 



T 

,2 / jj.j-2/ 



dtf'it). 



(46) 



In Appendix 1^ we rederive these expressions directly by 
second-order perturbation theory. 

A comparison between Eqs. (j45|) and (j46| reveals that 
in the weak- field limit, \ft\ ^ |A|, the ratio of the geo- 
metric and dynamical phases is equal to 



lAI' 



and thus independent of the field envelope. 

We conclude by estimating this ratio for typical exper- 
imental values. For instance, for the ls5(J = 2) — > 
2p3(>/ — 3) transition in argon [l^, the wave length 
A = 812 nm, the resonance detuning A = 3 • 10^ s~^, 
the velocity Vy = 700 m/s and angle a = 10~^, we obtain 
\"f / P\ = 0.1, which is feasible in an experiment. 



I3e 




cosh 9 



and according to Appendix [B] we find 



lAIr 



aarctan(a) — — ln(l -|- a^) 



(47) which for ^ 1 reduces to 



Z UJfy 



For the Eckart envelope the expression Eq. (|43| for 
the dynamical phase takes the form 



(51) 



(52) 



(53) 



In the last step we have recalled Eq. ([50)) for -fE- 

Hence, we confirm the fact that the ratio of the geo- 
metric to the dynamical phase is given by Eq. (j47p . 



VII. CANCELLATION OF DYNAMICAL PHASE 



VI. APPLICATION TO ECKART ENVELOPE 

In this section we evaluate the geometric and dynam- 
ical phases for the Eckart envelope defined by Eq. ([55]). 
For the details of the integrations we refer to Appendix 
mandO 



A. Geometric phase 

For the Eckart envelope we have T = oo and Eq. ([M)) 
takes the form 



7-E = WaT / d9 



cosh 6* 



V cosh^ 9 + a'^ 



(48) 



where we have introduced 9 = t/r and used the symme- 
try of the integrand. 

In Appendix [B] we perform this integral and find the 
geometric phase 



7b (x) = ~^UaT\n{l + a'^{x)) , 
which in the weak-field limit a{x) <C 1 reduces to 



(49) 



(50) 



We conclude by noting that we can rederive the ex- 
pression Eq. p9)) for by decomposing the path in the 
parameter space into a sequence of closed circuits and 
calculating the sum of the fluxes through each of these 
areas, as shown in Appendix [Cl 



Next we use the results obtained in the previous sec- 
tions to propose a scheme to cancel the dynamical phase, 
which always dominates the geometric one. Moreover, 
due to its dependence on the energy of the system, the 
dynamical phase is particularly sensitive to the slightest 
change of the parameters. 

From the expressions Eqs. (P7|) and ([50]) for the total 
phases acquired by the ground and excited states, we re- 
call that the dynamical part given by /3 sign( A) depends 
on the sign of the detuning A, whereas the geometric 
part 7 only on its absolute value. 

This fact allows us to suggest a rather intuitive scheme 
to compensate the dynamical phase. We propose to use 
two consecutive interactions of the atom with the stand- 
ing waves, that is, firstly with blue-detuned waves (A > 
0), secondly with red-detuned waves (A' = —A < 0). 
Here we have introduced a prime to indicate the second 
standing wave. As a result of the opposite signs of the 
detunings, the dynamical contributions cancel each other 
provided the condition k cos a — k' cos a' is fulfilled, or 
when the position-dependent Rabi frequencies n{x) and 
il'{x) defined by Eq. (fH]) are equal. 

At the same time, the geometric phases add up and 
result in the total phase 



<^f*)=(Pg + (^; = 7 + 7' 



(54) 



of the ground state. 

Thus, in a scheme of two consecutive scatterings of the 
atom by oppositely detuned standing waves, we obtain 
a cancellation of the dynamical phase and summation of 
the Berry phase. Of course, the time interval between 
the two interaction zones should be larger than the in- 
teraction time T itself, in order to be consistent with the 
adiabatic approximation on the turn-on and turn-off of 
the interactions. 
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VIII. FOCUSING DUE TO GEOMETRIC 
PHASE 

In the previous section we have shown that in our 
scattering setup with first an interaction with the red- 
detuned wave, and then with the blue-detuned wave, the 
atom acquires only the geometric phase. We now demon- 
strate that this phase manifests itself in a focusing of the 
atomic wave packet. 

For this purpose we assume for the wave function of 
the center-of-mass motion of an atom in the ground state 
moving in the y-direction a Gaussian 



vI/o(x) = 



1 



(V^Axo)i/2 



exp 



2Aa 



(55) 



of width Axo and the additional scattering-induced geo- 
metric phase jix). 

The time evolution of this wave packet in the absence 
of an external field is given by 



M 



CO 

j dx' exp 



'.M 
2hV 



„'^2 



*o(a;'). 
(56) 

In the case of the Eckart envelope and in the limit 
of AlAxo < 1 and (|rio|/|A|)(A;Aa;o) < 1, the geomet- 
ric phase 7(x) given by Eq. (j49|) near the nodes of the 
standing wave, for instance near x = 0, is a quadratic 
function 



h X 
7(a;) = -7: 



of X with 



2 Axl 



i^nT ■vMr(fc COS aAxo)"^. 



(57) 



(58) 



|AP 

Here we have used the definitions Eqs. (|14p and (|40| . 

By substituting the initial wave function Eq. (|55|) into 
Eq. (|56p . we obtain the time-dependent distribution 



\^{x,t)\'^ 



1 



V^Aa;(t) V Aa;2(t), 
of finding an atom with the coordinate x, where the width 



(59) 



Ax{t) EE Axo\ 1 - & 



(60) 



is determined by the initial width Aa;o and the Berry 
phase contribution. Here tg = MAxq /His the character- 
istic time of field-free spreading, that is 



Ax^'\t) = AxoJl + ^^. 



The minimal possible width 



AXri 



Axq 



(61) 



(62) 



A Ax(t) 








: 












Ax"" 


>- 



Figure 3: Influence of the geometric phase on the free prop- 
agation of an atomic wave packet represented by its time- 
dependent width Ax{t) given by Eq. (|60p . For an atom in 
the ground state (Asg and 6 = 5) the wave packet first fo- 
cuses and then spreads, whereas for the excited state {Axe 
and b — —5) it spreads from the beginning. This spreading is 
larger than that associated with the free propagation of the 
Gaussian (Ax'"' and b = 0) given by Eq. ((6T|). 



of the packet is reached at the time 

b 



1 + 62 



t. 



(63) 



In Fig. 3, we present the focusing of the atomic 
wave packet induced by the geometric phase for rela- 
tively small values of the parameter b. For example, the 
value 6 ~ 5 is achieved for WqT = Att, the Rabi frequency 
rio = 1.8 |A| and the initial width kAxo = 0.25. For 
b — 5 the expression Eq. (|62p for the minimal width 
predicts a shrinking up to around five times compared 
to the initial width, as indicated by the solid curve on 
Fig. 3 denoted by Aa;^. After the point t « 0.19 is the 
wave packet expands faster than the free wave packet 
represented by the dashed line Ax'-^-' . 

The atom in the excited state acquires the same geo- 
metric phase as the ground state but with the opposite 
sign. Therefore, the parameter b determining the nar- 
rowing appears with the opposite sign and results in an 
accelerated spreading rather than a focusing of the wave 
packet, as shown by Ax^ in Fig. 3. 



IX. CONCLUSIONS 

In this article we propose a scheme to observe the ge- 
ometric phase in atom optics based on the scattering of 
a two-level atom by two consecutive standing light waves 
with the same envelope but opposite detunings. The dy- 
namical and the geometric phases acquired by the two- 
level atom during the interaction are calculated within (i) 
the rotating wave approximation, (ii) the Raman-Nath 
approximation, and (iii) adiabatically slow switch-on and 
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-off of the interaction. Both the dynamical and the geo- 
metric phases are evaluated for a field envelope given by 
the Eckart function. 

We now specify the conditions under which these ap- 
proximations are valid and self-consistent. Indeed, the 
velocity acquired due to a resonant atom-field interac- 
tion can be estimated as 

where flo is a characteristic value of the Rabi frequency. 
Hence, we can omit the kinetic energy operator in the 
Schrodinger equation when {vx)t <^ or 

UJrecT'^^O ^ 1, 

with ojrec = hk^/{2M) being the recoil frequency. 

Moreover, we can neglect the spontaneous emission due 
to the small interaction time r, provided 

Ttw < 1, 

where F and w are the spontaneous emission rate and the 
occupation probability of the excited level, respectively. 
In the case of near resonance, |A| > fio, the maximum 
value of the population probability is ~ (flo/A)^. 

In a sequence of two such scattering arrangements with 
opposite signs of their detunings the dynamical phases 
compensate each other, whereas the geometric phases 
add. Therefore, the final state of the atom is different 
from the initial one only by the acquired Berry phase, 
which in our case depends not only on the internal, but 
also on the external degrees of freedom, such as the posi- 
tion of atom. This dynamical phase cancellation provides 
us with a completely different technique of measuring the 
geometric phase, employing the narrowing of the atomic 
wave-packet prepared initially in the ground state. This 
novel suggestion of the Berry phase measurement based 
on self-lensing beneficially differs from the previous inter- 
ferometric schemes of observing the geometric phase and 
might also constitute a useful tool in atomic lithography. 

The familiar WKB-technique is used as an independent 
method to verify the results obtained within the standard 
approach [l| and the results are shown to coincide. How- 
ever, we emphasize that the treatment of Ref. [l[ is more 
general in the sense that it can be employed for any de- 
pendence of the interaction envelope and the detuning on 
time. We consider in this paper the particular case of a 
constant detuning. 

We conclude by emphasizing again that the Raman- 
Nath regime takes place for a small interaction time t, 
that is r <C l/{u}rec^oy^^ ■ In this case, we are allowed to 
neglect the kinetic energy operator. However, for large 
interaction times, that is for l/(wreci^o)^^^ ^ r ^ l/F 
it should be taken into account. To do this in an exact 
way we could consider the scattering of the atom by a se- 
quence of two near-resonant running rather than stand- 
ing light waves [l^. In this case the formulae for the 
geometric and dynamical phases are analogous to those 



of the standing wave case, but independent on the exter- 
nal atomic degrees of freedom, i.e. the atomic position. 
Therefore, no lensing occurs and we can only use the 
standard interferometric scheme for the observation of 
the geometric phase. 
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Appendix A: Alternative derivation of phases 

In this appendix we present alternative routes toward 
the expressions Eqs. (p9l) and p3|) obtained within 
Berry's formalism. To solve the Schrodinger equation 
we make the ansatz 

|*> = Ae{t; r) e-'("-=-'^/2)*|e) + Ag{t; r) e-'("<'+'^/2)*|5), 

(Al) 

which is different from Eq. (O due to the effective de- 
tuning 

A = A + LOa- (A2) 

Substituting Eq. (jAll) into Eq. (U) and making the 
rotating wave approximation, we arrive at 

'4(t)^Kf-lA)(t)' 

where in contrast to V given by Ec^. (|12l) the coupling 
matrix element 

V(t; r) = Ml{x)f{t) = hn{x, t) (A4) 

is now real and depends on time only through the enve- 
lope /(i), describing the adiabatic turn-on and turn-off 
of the interaction. 

Two approaches to solve Eq. (jA3p and obtain the rel- 
evant phases offer themselves : (i) the WKB-technique 
and (ii) second-order perturbation theory. 

1. WKB-approach 

In order to transform the two first order differential 
equations, Eq. (jA3[) . into a single one of second order, 
we introduce the two functions 
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resulting in the differential equation [43 



1 / X 2 



n^ + 2i—]u^0 (A6) 
at I 



for u = u{t). 

This equation can be represented in a form similar to 
the stationary Schrodinger equation in position space by 
introducing the dimensionless variable ^ = i/r, where 
r is the characteristic time scale of the envelope /. In 
terms of ^, Eq. (|A6I) reads 



1^ 



~7 i dtl\ 

£ +7: — tt]u = 0, 



where 



(A7) 



(A8) 



Equation (jA7l) is analogous to the stationary one- 
dimensional Schrodinger equation where ^ and e play 
the role of the "effective coordinate" and "effective mo- 
mentum" , respectively. The small parameter 1 /t mimics 
Planck's constant in the conventional WKB-approach. 

Due to the adiabaticity condition 



1 de 



(A9) 



we can employ the semiclassical approach to search for 
the solution u(^) in the form 



(AlO) 



where the complex- valued function S{^) obeys the equa- 
tion 



dS_ 



T d^^ 



i dfl 
27 ~ 



(All) 



Within the WKB approach S{£,) is expanded into the 
perturbation series 



(A12) 



in powers of 1/t and from Eq. (|A11[) we find the zero- 
order term 



5(0)(0 = ± J d^eie) (A13) 



-T/r 



and the first adiabatic correction 
with 



lAI 



1 



A2 



1/2 



(A14) 



(A15) 



The second term on the right-hand side of Eq. (IA14P 
is a consequence of the last term in the left-hand side of 
Eq. (jAlip and appears in addition to the conventional 
WKB solution. Moreover, the next-order correction to S 
gives negligible contribution to u. 

The general solution of Eq. (|A7p can then be written 
in the form 




±iT J dem')T 



-T/t 



.(A16) 

where the coefficients A± can be found from the initial 
conditions Ag{-T) = 1 and Ae{-T) = for Ag{t) and 
Ae{t). By using Eq. (jASp and the connection between 
u and V we get u{£,) — 1/2 and {du/d£) — iKt/A at 
e — —T/t^ which results in 

A± = i(l±signA) = ie(±A). (A17) 

Here we have used the fact that at ^ = ^T/t the Rabi 
frequency 17 defined by Eq. (jA4p vanishes and therefore 
J\[{-T/t) = 1 and 2e{-T/T) = |A|. 

The same procedure can be applied to find v{t), leading 
to v{t) = —u{t). We then obtain the relations Ag = 
u — V = 2u and A^. — u + v = 0. The latter confirms 
the fact that there are no transitions into the excited 
state due to the adiabatically slow time dependence of 
the interaction. 

Thus, the probability amplitude Ag is given by 



±iT / dt'e{t')T 



InAf 



-T 



(A18) 

By taking into account the exponential prefactor 
exp(-iAV2) in Eq. we obtain from Eq. (IA18|) 

the total phase 



= / dt { V A2 + n^t) - |A| ) (A19) 



acquired by the ground state at t — T, when the interac- 
tion switches off, that is Ag{T) = exp(ii^g). 

According to Eq. (|A19p the total phase tpg is a function 
of the effective detuning A = A -|- Wq. When Wq, <^ |A| 
we can expand (pg into a Taylor series 

(^g(A) =(^g(A+t^„) =(^g(A)+W„^ (A20) 

over Wq/|A| and arrive at 



A 
"2" 



dt 



[Vl + a^P{t) 
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-T 



- 1 



(A21) 



Here we have recalled the definitions Eqs. (|40p and (jA4p . 

A comparison between Eqs. ([27|) and (|A2ip reveals 
that the first term in Eq. (jA21l) gives the expression Eq. 

for the dynamical phase, whereas the second term 
is the geometric phase given by Eq. (|39|) . Moreover, Eq. 
(|A20p shows that the Kramers-Kronig-like relation, Eq. 
(|44l) . between the dynamical and geometric phases is a 
consequaence of a Taylor expension. 

We conclude by noting that E g. (|A6[) for the Eckart 
envelope can be solved exactly [43| in terms of hyper- 
geometric functions. In the adiabatic limit, |A|r ^ 1 
and LUa |A|, the probability amplitudes Ag and 
obtained from the exact solution coincide [42| with those 
derived within the WKB approach. 



2. Perturbation theory 



The Schrodinger equation (|A3|) can be solved pertur- 
batively using the coupling matrix element V as the 
expansion parameter. Indeed, the second-order correc- 
tion to the probability amplitudes = A^ e*^*/^ and 
a„ = Ag e^'^*/^ are given by 



Similarly, we obtain 



dg{T)^ag{-T) + ~a'f'> ^~ag{~T) 



dt 



~ag{T)^~ag{-T)e'^. (A26) 

By using the definitions Eqs. ([40]), ((X2l) and (fM)) we 
derive for the total phase given by Eq. (|A25|) the ex- 
pression 



4 A-hWa 



dtf{t). 



(A27) 



In the case when cj^ ^ |A| the phase reads 

T T 

^ = ^A I dtfit) - I dtfit). (A28) 



With the help of Eq . (|27| we find that the first and 
second terms in Eq. (jA28[) give the dynamical and geo- 
metric phases in the weak-field limit defined by Eqs. p5)) 
and P5|) . respectively. 



Appendix B: Evaluation of integrals 



and 



dtV{t)e''^' / dt'Vlt') 



(A22) 



In this appendix we evaluate the integrals Eqs. 
and (|51l) determining the geometric and dynamical phase 
in the case of the Eckart envelope. We note that some of 
the relations are also useful for Appendix [Cl 



dtV{t)e''^' I dt'V{t')e 



(A23) 

By using the adiabaticity condition Eq. (fT9)) . we find 



dt'V{t')e^ 
which results in 

ae(T) = ae(-T) -1-5(2) =ae(-T) 



A 



1-1 [d& 

4 7 A 

-T 



that is 



with the phase 



de(T) = de{-T)e 



1 f , V^{t) 
4. J A 

-T 



(A24) 



(A25) 



1. Integral determining the geometric phase 

In order to calculate the integral 



Ie= I dO 





cosh 9 



_ vcosl? 
we recall the integral relation 
dO 



1 



(Bl) 



cosh 9 

V cosh^ e -f 

and find 

Ie = lim < In 



= = In (sinh6' + \/ cosh^ 6* + 
2 V 



sinh 9 + \J cosh^ 9 + 0? 



which with the help of the asymptotic behavior 

1 a 



in the limit of i 



sinh 9 ^ cosh i 



-> oo reduces to 
/£ = -^ln(l + a'). 



(B2) 
(B3) 

(B4) 
(B5) 
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2. Integral determining the dynamical phase 

The dynamical phase is determined by the integral 




1 



(B6) 



lE = a' 



1 



cosh^ e 1 + \/l + a2 cosh"^ 6 ' 

u 

The change of the integration variable 

/ a2 

sin I? = \ tanh^ 

V l + a2 

gives 

arctan(a) 



Ie = a 



arctan(a) — 



1 + Vl + a^ cos-!? 
When we recall the integral relation 

d^ 



A + Bcos'd 



In 



^/W^~W tan(i?o/2) + A + B 



yfW^~A? tan(i9o/2) - A-B 



for A < B, we arrive at 



Ie = aarctan(a) — - ln(l + a^) 



(B8) 



Appendix C: Evaluation of flux 

In this appendix we evaluate the flux through the area 
in parameter space defined by the trajectory following 
from the Eckart envelope of the electromagnetic field. In 
this case the separation of the trajectory from the origin, 
determined by the strength of the envelope, first increases 
from zero to a maximum value and then decreases again 
to zero. As a result, we obtain infinitely many windings 
and thus infinitely many areas of different sizes. 

In order to calculate the total flux through them, we 
present the expression Eq. p7p for the geometric phase 
7£; as a sum 



El 



pi4>) 



pdp 



(1 + p2)3/2 



E7 



(m) 



(CI) 






Y 


/ 

1 


\ 

\ 


• ^ ni 

(B7) 









Figure 4: Translation of the field envelope / = f{t) during 
two symmetrically located time intervals (top) into a closed 
circuit in parameter space (bottom). Indeed, the m-th circuit 
arises from the envelope / = f{t) during the time intervals 

and t,n-i < t < tm- 



(C2) 



of the fluxes through the m-th area with 

COSh[(/)/ (uJaT)] 

Here we have recalled the definitions Eqs. (p4|) 
and PI|) . 

In order to perform the integration in Eq. (ICip we 
need to determine the angles corresponding to the path 
defining the m-th area. This path is dictated by the 
envelope during the time intervals —tm < i < —tm-i 
and tm-i < t < tm, where tm = rmi/uja, as shown in 
Fig. m Since = ujat these time domains translate 
into the integration intervals —4>m ^ 4' ^ —4'm-i and 
0m-i < (l> < 4>m with 0m = mn. 

Therefore, the geometric phase ^^"^"f given by the flux 
through the m-th area reads 



-(m-l)7r p(</)) 



pdp 

(l+p2)3/2 
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rriTT p(0) 

1 f . . f pdp 



2 J ""7 (l+p2)3/2 ' 

{m~l)-!T 



(C3) 



which due to the symmetry of the Eckart envelope giving 
rise via Eq. (jC2p to p(— 0) = p{(j)) simpHfies to 



v(™) 



d(j) 



(m-l)7r 



pdp 



(1 + /52)3/2 ' 



(C4) 



or 



7^-^=-.r / dol ^ -1 (C5) 

V V cosh^ 9 + I 



with = rmr / (lOaT). 

With the help of the integral relation Eq. (jB2[) we 
arrive at 



(™) = War In 



where 



7 



i^m = sinh( 



Frn TT 



cosh 9m + a? 



(C6) 



through all areas, that is 

°° ^ f F TT 

= 7^™'' = lim I In — — 



Using the functional relation 



(C7) 



In- 



F„ 



\nFm - lnFm_i 



of the logarithm we find 



7s=Warlim (InFAr - InFo - TV— ) . (C8) 

N^oo \ Ll!„T ' 



The asymptotic expansion Eq. (|B4|) yields 



^ exp ( N— ) , (C9) 



which together with Fq = (l + a^Y^^ leads us to the total 
geometrical phase 



IE = -z- m ( 1 



(CIO) 



The total geometric phase is the sum of the fluxes which coincides with Eq. (|49 
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